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INTERACTIVE COMPUTER MODELING OF COMBUSTION CHEMISTRY 

and 

COALESCENCE-DISPERSION MODELING OF TURBULENT COMBUSTION 

David T. Pratt 
University of Washington 


The goals of this research project are as follows: 

1. Develop an interactive computer code for simulation of a high- intensity 
turbulent combustor as a "single point" inhomogeneous stirred reactor [l], This 
will be developed from an existing batch processing computer code CDPSR [2]. 

2. Use the interactive CDPSR code as a guide for interpretation and direction 
of DOE-sponsored companion experiments utilizing Xenon tracer with optical laser 
diagnostic techniques to experimentally determine the appopriate mixing frequency, 
and for validation of CDPSR as a mixing-chemistry model for a laboratory jet- 
stirred reactor. 

3. Incorporate the coalescence-dispersion model for finite rate mixing into an 
existing interactive code AVCO-MARK I, to enable simulation of a combustor as a 
modular array of stirred flow and plug flow elements, each having a prescribed 
finite mixing frequency, or axial distribution of mixing frequency, as appropriate. 

4. Further increase the speed and reliability of the batch kinetics integrator 
code CREKID [3] by rewriting in vectorized form for execution on a vector or parallel 
processor, and by incorporating numerical techniques which enhance execution speed 
by permitting specification of a very low accuracy tolerance [4]. 
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*** HARK2I *** 

•HARK2I* IS AN INTERACTIVE VERSION OF THE MARK-II COMBUSTOR MODEL. 
THIS IS A PRELIMINARY DESIGN TOOL, A MEANS OF GAINING INTUITIVE INSIGHT 
INTO EFFECTS OF CHANGES IN FUEL-AIR MIXING OR PARTITIONING ON TURN-DOWN 
RATIO, COMBUSTION EFFICIENCY AND POLLUTANT FORMATION RATES. 

AN INITIAL DATA SET IS TAKEN FROM DATA FILE * HAKK2.DAT* BUT CAN BE ALTERED 
INTERACTIVELY, AND USED IN CONSECUTIVE RUNS. 

HAKK-II REPRESENTS A SIMPLE BRAGG COMBUSTOR CONSISTING OF A 
MAXIMUM OF 9 FLOW ELEMENTS WITH THE ADDITION OF A SINGLE RECTCLE 
ELEMENT. FLOW ELEMENT TYPES MAY INCLUDE: 

1) NON-REACTING MIXERS (‘MIX’), IN WHICH THE CHEMICAL REACTIONS ARE 
ASSUMED TO HAVE STOPPED DURING THE MIXING PROCESS, 

2) PERFECTLY STIRRED REACTORS CPSR*), WITHIN WHICH INTENSE SELF- OR 
BACK-MIXING IS ASSUMED TO OCCUR, SO THAT THERE ARE NO AXIAL GRADIENTS! 

3) PLUG FLOW REACTORS < *F‘FR* ) . 

THE USER MAY DEFINE THE MODEL AS HAVING UP TO 9 ELEMENTS IN 
SERIES WITH AIR AND FUEL INLET JETS AT EACH ELEMENT. THE RECYCLE ELEMENT 
MAY BE OF ANY OF THE THREE FLOW TYPES, AND MUST RECYCLE FROM A HIGHER 
NUMBERED ELEMENT TO A LOWER. COOLING BOUNDARY LAYER EFFECTS AND CHEMICAi 
REACTIONS WITHIN THE BOUNDARY LAYER ARE NOT CONSIDERED. 

PLEASE WAIT A MOMENT WHILE INITIALIZATION IS COMPLETED. 

INITIALIZED — PRESS ''RETURN', TO BEGIN — 
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ELEMENT ) 
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RECYCLE 


(SO. IN) 
1.A400E+02 
1.4600EK52 
1.4A00EV02 
l.ASOOEtOO 
1 . 4A00F ►02 


LENGTH 

FLOW 

< INCHES) 

TYPE 

4.0000E-01 


1 • OOOOt-O l 

h LX 

1 » 5000E+00 

PSR 

2 . OOOOE-0 1 

MIX 

1 .OOOOE+OO 

MIX 


INLET AIR 
<i.:m/s> 
1.3-.rOELOO 
A.9SA0E-01 
O.OOOOE-Ol 
7.VLQ0E-O1 
RECYCLE 20 
OUTFLOW TO 


INLET FUEL 
( LDM/S > 
9.7200E-02 
O.COOOE-O l 
O.OOOOE-Ol 
O.OOOOE-Ol 
OOX Or 5 3 
i 2 INFIOU 


AIR TFMP = 2 . 1 OOOt 102 F 
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COMBOS TOR PRESSURE • 2.?100E A 00 AIM 

LOWER HEATING VALUE •> 1 . RSOOErOA STU/LBn 


s: l 2CT AN OPTION BY NUMBER: 

0- RUN WITH THIS DATA SET 
-X - CHANGE AIR TEMPERATURE 
-? • CHAN- E FUEL TEMPERATURE 
-3- CliAr-uE LOWER HP' AT I MO VALUE 
rOTIO.-P (0-7) 7 


-A- CHANGE NOMINAL CC.1BUST0R PRESSURE 
-5- CHANCE RECYCLE CLEMEN f STATUS 
-6- C!'-NC..' FLO'J CLThENTS STATUS 
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